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Abstract 

In models of universal extra dimensions, gravity and all standard model fields propagate in 
the extra dimensions. Previous studies of such models have concentrated on the Kaluza-Klein 
(KK) partners of standard model particles. Here we determine the properties of the KK gravitons 
and explore their cosmological implications. We find the lifetimes of decays to KK gravitons, 
of relevance for the viability of KK gravitons as dark matter. We then discuss the primordial 
production of KK gravitons after reheating. The existence of a tower of KK graviton states makes 
such production extremely efficient: for reheat temperature Trh and d extra dimensions, the energy 
density stored in gravitons scales as T^^^^"^ . Overclosure and Big Bang nucleosynthesis therefore 
stringently constrain Trh in ah universal extra dimension scenarios. At the same time, there is a 
window of reheat temperatures low enough to avoid these constraints and high enough to generate 
the desired thermal relic density for KK WIMP and superWIMP dark matter. 

PACS numbers: ll.lO.Kk, 12.60.-i, 95.35. +d, 98.80. Cq 
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I. INTRODUCTION 



Since the initial work of Kaluza and Klein yj, models of extra dimensions have played 
an important part in particle theory. These ideas have been revived in models of TeV-scale 
universal extra dimensions (UED) j^l^l^. In these models, as in Kaluza-Klein (KK) theory, 
both gravity and standard model (SM) fields propagate in D = 4 + d dimensions. 

In models where the extra dimensions are tori, every known particle has a tower of 
KK partner particles, each carrying KK number. Momentum conservation in the extra 
dimensions implies KK number conservation. Such models predict a plethora of unseen 
massless particles and a tower of stable KK particles, both of which are phenomenologically 
problematic. 

For more general geometries, unwanted massless modes may be projected out, and there 
are fewer or no stable KK particles. A particularly interesting situation occurs when there is 
only a Z2 symmetry (called KK-parity) in the extra dimensions. An example of this situation 
is the case D = 5 where the extra dimension is the orbifold S^/Z2, where is a circle of 
radius R For SM fields, the geometry projects out half of the massless modes, leaving 
only the 4D SM degrees of freedom, and also breaks KK number conservation, rendering 
almost all KK partners unstable. However, the geometry preserves KK-parity, and so the 
lightest KK particle (LKP) is stable, making it a viable realization of dark matter with 
an extra-dimensional origin Among SM KK partners, the LKP is often B^, the KK 
partner of the hypercharge boson |^. For mxK ~ TeV, is an excellent weakly-interacting 
massive particle (WIMP) dark matter candidate, with a thermal relic density consistent 
with observations and promising prospects for detection B, Isj IE jl^ 111 j^L, 

These studies, as well as those exploring the collider ja lisl l and low- 

energy 3, H 0, H H H H 113 implications of UED, have focused on the SM KK 



spectrum, neglecting the gravitational sector. KK gravitons couple extremely weakly, so 
their effects are insignificant for colliders and low-energy experiments. Cosmologically, how- 
ever, KK gravitons may play an important role. For example, the LKP may be the lightest 



KK graviton [26|, |27[. It is not easy to verify this, since loop diagrams involving gravi- 
tons are divergent. On general grounds, however, one expects the mass corrections to be 
5m/ m ~ (Mpii?)~^, since gravitons decouple as Mpi 00. This is potentially very small, 
and, since loop corrections to other KK particle masses are typically positive it is not 
unlikely that in UED theories, the lightest KK graviton is the LKP. In this case, all SM KK 
particles eventually decay to the KK graviton, and the KK graviton is the only possible KK 
dark matter candidate. 

KK graviton dark matter interacts only gravitationally and is superweakly-interacting 
massive particle, or superWIMP, dark matter |26l . |27| . As such, it is impossible to de- 
tect in conventional dark matter experiments. However, its potential impact on Big Bang 
nucleosynthesis (BBN), the cosmic microwave background, and the diffuse photon flux pro- 
vides alternatives for dark matter searches. In evaluating such signals, it is important to 
have accurate, rather than just order of magnitude, results for decay times of WIMPs to 
superWIMPs. We determine these in this study. 

In addition, KK gravitons may be produced in the early universe during the reheating era 
following inflation. The phenomenon of gravitino production after reheating in supersym- 
metric scenarios is well-studied HQ. The case of UED is qualitatively different, however. 



as there is an infinite tower of new particle states that may be populated at high tempera- 
tures, with the density of states growing rapidly at large masses, especially for large d. As we 
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will see, in UED primordial KK graviton production is in fact very efficient, and constraints 
on dark matter abundances and BBN provide stringent bounds on early universe cosmology 
in all UED scenarios, irrespective of which particle is the LKP and other spectrum details. 

The paper is organized as follows. In Sec. [H] we derive the interactions of KK gravitons 
coupled to UED SM fields in general toroidal dimensions. In Sec. IIIII we then obtain the 
KK graviton interactions for UED orbifold compactifications, specializing to the case of the 

/Z2 orbifold discussed above. For this case, we first show that boundary terms coming 
from the ends of the interval are typically negligible and then determine the spectrum and 
bulk interactions of the orbifold theory. In Sec. IIVI we find the decay rates for next-to- 
lightest KK particle (NLKP) gauge bosons 5^ ^ + 7 and fermions ip^ ^ + ip, of 
relevance to superWIMP dark matter. We then estimate the primordial abundance for the 
tower of KK gravitons produced after reheating in Sec. |3 and derive constraints on the 
reheat temperature in Sec. IVll Our conclusions are given in Sec. IVllI 



II. GRAVITON INTERACTIONS FOR TORUS COMPACTIFICATIONS 



To discuss graviton cosmology, we must first derive the interactions of KK gravitons in 
UED. Graviton interactions in extra dimensions have been discussed previously, particularly 
in the context of theories where SM fields are confined to four spacetime dimensions. (See, 
e.g., Refs. js^HH; we will follow the analysis of Ref. in discussing the linearized grav- 
itational action.) Here we instead couple the gravitational sector to extra-dimensional SM 
fields and also implement the orbifold projection. 

We begin in this section by determining the couplings of KK gravitons and KK SM fields 
in the case of torus compactifications. For D spacetime dimensions, the action for matter 
coupled to gravity is 



S 



D 



ME-' 



d'^xJlgl R + C{(j)) 



(1) 



where M^i is the Z)-dimensional Planck scale and C is the Lagrangian for all matter fields, 
which we denote generically by 0. We use coordinates x*'^ = (x'^, y*), where < < 27cRi = 

u. 

Linearizing about fiat space using Qmn = Vmn + hMN, where 77 = (1,— 1,-1,...,— 1), we 
find the Fierz-Pauli action 13^ 



oD 
'-'lin 



D-2 
D 



K{h)+Co{4>) + \hMNT''''{4>) 



(2) 



where the graviton kinetic terms are contained in 



K{h) ^- {h 



MN,PT T,MT c)TMN T ,Pj_oh h^'^' 



N 



with h = h^M and /i*^^ = r]^^^7]^^hpQ. Also 
is the fiat-space matter Lagrangian, and 



MN , 



dC 



'dg 



(0) 



MN , 



(3) 



(4) 



(5) 



9=V 
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is the flat-space stress-energy tensor. To reproduce the standard gravitational interactions 
among massless modes, we require Hi Li = {16ttGn)~^ = M|. Defining h - M^h and 

(p = M40, we find the hnearized action 



2M, 



(6) 



where the kinetic terms are now canonicaUy normahzed. 
To write this in terms of 4D fields, we define 



huN — 
Tmn — 



T T ■ 
T T 

-'-IV -i-ij 



(7) 
(8) 



where = For simplicity, we will henceforth take Ri = R for all i. Each component 

field is then conventionally decomposed into 4D fields through the Fourier expansion 



n 
n 



(9) 



and similarly for the stress-energy tensor component fields. Here n ■ y = J2i i^iU^, where 
n = {rii, . . . , Ud) and each rii runs from — oo to oo. We will also use the variable 



rrir, 



i?2 



1/2 



(10) 



At each massive level n, part of the fields A^j, A^^, and 0^- are absorbed to give mass to the 
graviton state h^^ in an extra-dimensional realization of spontaneous symmetry breaking. 
The resulting physical states are (defining ff = d'^^'Ukrii) 



h 



An 



/ t 



- S 



o ' J 
-2 



28^ 
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where we have defined 



and a satisfies 3{d — 1)0^ + 6a 



P 



_ njUj 

ij -»9 



(11) 

(12) 
(13) 

(14) 
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The Lagrangian for the physical fields needs to be treated slightly differently for the zero 
modes versus the non-zero modes. For the zero modes, the action is 



lint 



where 



'"hiat 



i i i 



2M 



rpOij 



The non-zero modes have the kinetic terms 

Sli^-^- = J d'x [Cl + Ci + L% 



where 



rf. 



'"A 



1 nr^....r ^ 1 or^r ^ 



and F^^j = d^A^j^ — d^A^^. The interaction terms of the 



on int 
'-'lin 



non-zero modes are 



where 



rn 

'"hint — 



pn p 'Xn{i,jrp—ni 



nn 



M4 

1 ~ 



±nrp-nij 



with 



rpn 



P 



(15) 



(16) 
(17) 



(18) 



V2(l + a(d- 1)) 



- V2P,,PjiT^i 



(19) 
(20) 
(21) 

(22) 

(23) 
(24) 
(25) 

(26) 



Note that in all expressions above, both (ni, n2, ■ ■ ■ , Ud) and (— ni, — n2, . . . , —rid) arc to be 
included in sums over n 0, and all (P pairs are to be included in sums over ij. 
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III. GRAVITON INTERACTIONS FOR ORBIFOLD COMPACTIFICATIONS 



The analysis so far is valid for compactification on a torus S'^. For reasons described in 
Sec. m we are most interested in the case where KK number is broken to KK-parity. This 
may be achieved by compactifications on orbifolds. 

For simplicity, we now concentrate on the D = 5 case with coordinates = {x^,y), 
where the extra dimension is an interval /Z2 of length vri? j^]- This geometry can be 
described as an orbifold of a circle compactification, where we orbifold by the symmetry 
y — > —y. In the SM sector, this symmetry is accompanied by the transformations — 
and V5 — V5 for the gauge fields. The 4D scalar is projected out, and the massless sector 
includes only the 4D gauge field. A similar projection on fermions removes half of the degrees 
of freedom, leaving only the chiral fermions of the SM. 

In the gravitational sector, we will project by the action h^^ — > h^^, h^^ — > — /i^s, and 
^55 ^55 under y — ^ —y. At the massless level, this preserves the 4D graviton 
while removing the gravi- vectors h^^^{x) and h!l^{x). The gravi-scalar h!l^{x) is preserved 
by this projection; we assume that some other physics stabilizes this mode and generates a 
large mass for it. (Th e phenomenology and cosmology of an extremely light gravi-scalar is 
discussed in Refs. |33|.) The final massless sector is, then, just the SM plus the 4D graviton. 



A. Boundary Interactions 

Prior to discussing the bulk interactions, we need to consider the boundary terms. These 
arise because the space we are considering is singular at = 0, vri?, and hence we expect 
new interactions at these points. The nature of these interactions can only be derived from 
an underlying microscopic theory (such as string theory) that smooths out the singularity. 

Nevertheless, we can make some qualitative statements about these terms. If the singular 
points are smoothed out to a size Ip-, these interactions are localized near the singular points 
in a region of size Ip- We can formally write these interactions in the form 

^-/"'^r^Ka^/i^)]-'^)' 

where f{w) goes to zero for w 1. We have assumed that the resolution of the singularity 
is such that KK-parity is preserved. 

The basic assumption we will make is that Ip -C vri?. This is not unreasonable; for 
example, a natural guess for the scale Ip would be the (higher-dimensional) Planck scale. If 
we then consider R~^^ ~ 1 TeV, we find ~ 10^° TeV, so indeed Ip <^ vri? in this case. If 
Ip ^ 7ri?, we can estimate the boundary terms as 

5~-^/(0) Id'xCix) , (28) 
Tin, J 

that is, the boundary terms are suppressed by a factor Ip/nR, which is small. Thus, to 
leading order, we can ignore these terms. 

KK number- violating terms are generated at one-loop [s^l- The tree- level spectrum 
already breaks KK number, since certain modes are projected out by the orbifold action. 
These are translated at loop-level to interactions that violate KK number. The typical decay 
widths of the higher KK modes are therefore ~ arrin- 
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There is another source of boundary terms: certain manipulations of the bulk terms 
require integration by parts, which produce boundary contributions. These terms are also 
suppressed for the reason given above. 



B. Bulk Interactions 



Because we project out modes that are even or odd under y —y, it is more convenient 
to replace the conventional Fourier expansion of Eq. © by the expansion 



n>0 
n>0 

T,,ix,y) = Tl,{x) + V2Y. 



n>0 



Kt^uix) cos — + h_^^{x) sm — 

"^+5.(^)cos^ + A:!,,(x)sin^' 

,„ , , ny ,„ / N . ny' 
0+ 55 (^) cos — + 0„ 55 (^) sm — 



n>0 



T+^,Ax) COS — + T_^^(x) sm — 



(29) 
(30) 
(31) 
(32) 



and similarly for the other component fields. All /i"^^, ^+5;^, ^+^57 and 01 55 fields are 
projected out by the orbifold. 

The physical graviton field at KK level n > is 



55 



In this new basis, the linearized action for the G fields is 



J] 



lin 



d X 



n>0 



where 



2M, 



Gint 



n rjiTL fll/ 



GIX 



2M: 



The KK graviton propagator is^ 



k"^ — mi + ie 



(33) 

(34) 

(35) 
(36) 

(37) 
(38) 

(39) 



^ This analysis differs by a factor of 2 from the journal version of Ref. [s^, but agrees with the later e-print 
version hep-ph/ 9 8 11 3 5 Oy 4 . 
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where 



k^kp 



mt 



mi 



' + 2 - 



ki/kp 



mt 



kpky 
ml 



kpkfj 
ml 



(40) 



and mn = n/R. 

To determine the KK graviton couphngs, we must determine the matter stress-energy 
tensor components t!^^'^. For gauge fields, 



Tmn — Fm^Fpn — -^VmnFpqF^^ . 
We expand the gauge field in harmonics 

a^(x,y) = al^(x) + V2j2\al^{x)cos'^ + a''_^{x)sm'^ 



n>0 



R 

ny 



R 

ny 



a5{x,y) = a%(a;) + V2^ a!;: 5(2;) cos — + a!^ 5(0;) sin — 



n>0 



R 



(41) 

(42) 
(43) 



All a"^ and a" 5 fields are projected out by the orbifold. 
The physical gauge field A'^ at KK level n > is 



R 



n 



and identifying coefficients of cosny/R, we find 



p,v 



E 

m=0 



F;-ppn-m _ ^^^^j,n^j,n-mpa ^ rn,,m„{A^ A^"^ - ^r^^^A^A^-^ ") 



(44) 



(45) 



where FJJl = 9^^^ - d^A"^. The G'^^(q)A'^(ki)A'^-'^(k2) vertex is therefore 

i r 

Xf,uai3 = TTVT 'napkip,k2v " rjtiakipk2u - r)i'i3kink2a + VmVi'lsi^i • ^2) 
ZlVl/^ L 

1 1 

-o'?/^;^ (^/^/^(/^l ■ /^2) - ki(ik2a) + mnmn-miVnaVufS " T^VnuVafs) + (a ^ /3) 



,(46) 



The fermion couphngs may be calculated similarly. The stress-energy tensor for 5D Dirac 
fermions is 



(47) 



where m^o is the 5-dimensional mass of the fermion. We will consider a Z2 action on 
the fermion which acts as ip{x,y) = —'y5ip{x,—y). This preserves the zero mode of the 



8 



left handed fermion, but not that of the right handed one. In this case, the expansion in 
harmonics is of the form 



ra>0 



n>0 

n 
L 



^l(^) cos — + ii)R[x) sm — 



, „ , , . ny . ,„ f \ 
ipAx) sm — + tipn[x) cos — 



R 



R 



The stress-energy tensor for fermions is 



E 

m=0 



(48) 
(49) 



(50) 



where ipfi{x) = 0. The KK graviton-fermion-fermion vertex for G'^^,{q)ip'^{ki)ip'^ ^{k2) 
interactions is then 



Y, 



AMa 



Note that if the fermions in this vertex are on-shell, only the last two terms remain. 



(51) 



IV. NLKP DECAY WIDTHS INTO LKP GRAVITONS 



Using the vertices of Sec. IIII B[ we can now calculate Feynman diagrams involving gravi- 
tons and SM fields. We are particularly interested in the decay widths for applications to 
models of dark matter [26, 27], where an NLKP SM field decays to a stable LKP graviton. 
The NLKP may either be a gauge boson, such as B^, or a fermion, such as r^. These decay 
widths may be calculated using the usual trace techniques or by helicity amplitude methods. 
We have done both and checked that they yield identical answers. We present the helicity 
amplitude analysis here. 

The polarization vectors of a massive graviton, given the normalization conventions of 
Sec. irrm are 



.±2 
.±1 



9 



Here and e° are the polarization vectors of a massive gauge boson; for a massive vector 
boson with p^^ = {E, 0, 0,p) and mass m, 



4(p) = ^(0,1,^,0) , 
= ^(0,-1,2,0) 



m 



ip,0,0,-E) 



(52) 
(53) 
(54) 



The graviton polarization vectors satisfy the normalization and polarization sum conditions 

, (55) 



E er,„e 



B 



flu per 1 



(56) 



where B^^p^ is as given in Eq. (jlOI)- 

With these conventions, the helicity amplitudes for the decay of a KK hypercharge gauge 
boson to a KK graviton are 



M[B'Hk)^G^\p)+^Hq))\ = cosewX^^'^^e^{k)ep{q)el:{p) 



cos 6w m^B^ 



m 



M{B'\k)^G'^\p)+^^{q))\ = cose^X'^-^6°(fc)eJ*(g)e2*(p) 



cos 9w m^i 



1 - 



m 



M{B'^ik)^G'^\p)+^^{q) 



V2M4 rriGi 
cosdwX^^-^et^{k)ep{q)ef,\p) 



cos 



w 2 



1 - 



51 



(57) 



(58) 



(59) 



The total squared amplitude, averaged over the initial three polarizations and summed over 
final states, is 



\ Y: \M{B'\k)^G"^'ip)+^'"{q) 



h,h'h" 



cos^ 9w m%i 
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Ml 



1 ^ 



m. 



6 + 3- 



m 



+ 



m 



51 



m 



Gi 



Gi 



and the decay width of a KK hypercharge gauge boson to a KK graviton is therefore 



rmi^r^ioA cos^ m|;i 
1 (i^ ^ G 7) = ,,,_,^o_4 



1447rM|m^i l 



Gi 



3 r 



1 + 3^^ + 64^ 



(60) 



(61) 



The decay width of a chiral KK fermion to a KK graviton may be calculated similarly in 
the helicity amplitude formalism. For = {E, 0, 0,p), the helicity spinors are 



u (p) 



^E + mj- 



(62) 
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where = (1, 0) and ^ ^ = (0, 1), and we take the Dirac representation. 

The hehcity amphtudes for a chiral fermion decaying to a KK graviton are, then, 



M {fl%{k) 



Y^''uT{q)u^{k)ef^*{p) 



m 



/I 



2M4 rriQi 



m 



m 



3/2 



Y^^u^{q)u^{k)e%{p) 



m. 



m 



/I 



3/2 



(63) 



(64) 



where we have assumed a massless zero mode fermion. The squared amphtude, averaged 
over the initial two polarizations and summed over final states, is 



h,h'h" 



1 m. 



6 Ml 



1 - 



3 + 2- 



ra 



ra 



(65) 



Gi 



The decay width of a chiral fermion to a KK graviton is 



1 



967rM? 



Gi L 



Gi 



4 r 



2 + 3 



Gi 



2 



(66) 



The decay widths of Eqs. (jfjlj) and provide accurate expressions for NLKP lifetimes 
in the LKP graviton scenario. For NLKP and LKP masses at the weak scale Mweak, and 
assuming mass splittings of the same order, the naive expectation of lifetimes r = F^^ ~ 
47rM|/(Mweak)^ ~ 10^ s - 10^ s is born out. NLKPs therefore decay after BBN, and their 
decay products are subject to rather stringent BBN constraints. 

Even relatively mild degeneracies may disrupt these expectations, however, as for Am = 
"^NLKP — ""^LKP ^ "^LKP the lifetimes scale as r oc (Am)^'^ and r oc (Am)"*^ in the gauge 
boson and fermion cases, respectively. This behavior is mirrored in the analogous case 
in supersymmetr y o f superpartners decaying to gravitinos. The supersymmetric case was 
analyzed in Ref. j27|, where alternative signals of super WIMP dark matter were identified. 
The similarity of Eqs. and (jMj) to the supersymmetric case indicates that KK graviton 
dark matter is also a viable superWIMP candidate with promising possibilities for detection. 



V. REHEATING AND KK GRAVITON PRODUCTION 

KK gravitons are produced copiously after the Big Bang. Inflation dilutes these gravitons 
away, but their number density may be regenerated during reheating. The situation is 
analogous to the case of supersymmetry, where gravitinos may be produced after reheating, 
and constraints bound the reheat temperature to Trh ~ 10*^ — 10^° GeV jisi 2^|. As we 
will see, however, the presence of a tower of KK levels leads to much stronger production of 
gravitons, and much stronger bounds on reheat temperature in UED models. In this section 
we estimate the density of KK gravitons produced after reheating for UED models with an 
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arbitrary number of extra dimensions d. In the next section, we discuss the cosmological 
significance of these results. Throughout we assume that the extra dimensions remain fixed 
in size. 

UED models are characterized by two hierarchically separated scales: the KK mass scale 
^KK ~ TeV, and the 4D Planck scale M4 ~ 1.7 x 10^^ GeV. For this analysis, cosmology 
introduces another scale, the reheat temperature Trh. As UED models are 4D effective 
theories of some higher- dimensional theory, they are valid only up to some cutoff scale Mg, 
which is much smaller than M4 |35||]. We therefore assume Trh < Mg, and so Trh M4. 

In terms of these scales, the expansion rate of the universe is if ~ T'^/M^. The interaction 
rate of SM particles and their KK partners with each other is asun ~ a^T. The decay rate 
of SM particles at KK level n may also be estimated to be Fsm ~ arrin- Given the hierarchy 
between T and M4, 

(^sun , TsM > H , (67) 

and so SM particles and their KK partners remain in thermodynamic equilibrium as the 
universe cools after reheating. In contrast, the time scale for graviton interactions is very 
long. For example, the interaction rate for (f)a~^ 4'b ^ + where k, I, m, and n specify 
KK levels and a, b, and c label SM degrees of freedom, is acu ~ T^/M^. Similarly, as 
discussed in Sec. IIVI the decay rates of KK gravitons are Tq ~ ml^/M^. These rates are 
therefore far below the expansion rate, 

aGn,TG<t:H, (68) 

and gravitons never reach thermodynamic equilibrium. 

The qualitative picture, then, is that after reheating, the SM degrees of freedom exist in 
thermal equilibrium. Occasionally, they produce KK gravitons, which are meta-stable and 
decay long after Big Bang nucleosynthesis. If overproduced, they may overdose the universe 
or their eventual decay products may destroy the predictions of Big Bang nucleosynthesis. 

We now work toward a more quantitative estimate of graviton abundances. During the 
era of graviton production, the expansion rate H is given by 

H = -^Pn = Y^^9.{T)T , (69) 

where pn and g*{T) are the total energy density and the effective number of light degrees of 
freedom, respectively. The entropy density is 

5 = -^9*iT)T' , (70) 
and the number density of a massless bosonic degree of freedom is 

no = . (71) 

TT 

For UED theories, 

g.{T)=g^^D,{T) , (72) 

Here gf^ is the effective number of degrees of freedom per KK level, and is a model- 
dependent constant. In all models, however, there are more degrees of freedom at excited 
KK levels than at the zero mode level, where many degrees of freedom are projected out. For 
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the SM, the effective number of degrees of freedom, organized by spin, is gf^ = go^ + lgi^ + 
gSM ^ 4 iQQ ^ 24 = 106.75. For the S^/Z2 UED model, at each KK level > a gauge 
boson field has 3 degrees of freedom and a fermion field has 4 degrees of freedom, so the 
total effective number of degrees of freedom is g^^ = fl'o^^ + |5'i/2| = 197.5. A similar 

counting may be done for other UED models once the number of additional dimensions and 
the orbifold is specified. 

The function Dfi{T) counts the number of excited modes in the thermal bath at tem- 
perature T. Dct{T) may be approximated by counting all modes with mass below T. This 
yields, for d extra dimensions. 



rriKK 



(73) 



where 



Trd/2 4 1 

V,= ^^ =2,7^,|7^,^7r^... , (74) 



r(i + 



d] "'"'3 '2 



for d = 1, 2, 3, 4, . . ., is the volume of a unit spherical ball in d dimensions, and the factor of 
1/2"^ in Eq. (|7n|l accounts for the restriction to n with non-negative components. 

The number density of KK gravitons at level determined by the Boltzmann 

equation 

+ mncn = Con , (75) 



dt 
where 

is the collision operator. All graviton destruction processes are negligible, given the low 
abundances and decay rates of gravitons. We will parameterize the collision operator as 

CGn=Ca[g4T)nof , (77) 

where 

Here as is the strong coupling constant and C is a constant. If every light degree of freedom 
could interact with every other light degree of freedom with cross section a, Cqh with 
C = 1 would be a reasonable estimate. However, global and local symmetries restrict which 
reactions are possible, and not all interactions involve strongly-interacting SM particles. A 
detailed calculation of C can be done once a specific UED model is chosen and its spectrum 
determined. Here we keep the analysis general by leaving C as a free parameter. Based 
on the results of detailed studies of gravitino abundances from reheating in supersymmetric 
models [i^, we expect values of C ~ 0(0.01), and we will consider a range of 0.001 ^ C ^ 
0.1 in present numerical results below. Given the high power dependence of the graviton 
abundance on Trh, we will see that bounds on Trh are rather insensitive to even this generous 
range for C. Note also that, as C = 1 is certainly an overestimate, if a scenario passes all 
constraints even with C = 1, it is certainly allowed. 

With these definitions, the KK graviton number density satisfies 



dt 



+ ?>HnGr. = Ca [^*(T)no] . (79) 
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This is most conveniently solved by changing variables Hqu Yqr = Hqh / s and t T. 
Adiabaticity implies that entropy S = sR^ oc g^{T)T^E? oc T^^'^E? is conserved, and so 



Ids _ ^IdR _ dT 
s dt R dt ' dt 



3 + 



-HT . 



With the relations of Eq. ()80p. the Boltzmann equation becomes 

dYcn 3 + d 1 



dT 



3 HTs 



Ca [g^{T)noy 



(80) 



(81) 



Yqu changes until production stops at temperature T ~ m„ and then remains constant 
until gravitons begin to decay. After BBN and before KK gravitons decay, then. 



Trh 3 + d 1 2 
dT Ca [g,{T)no] 

nin J II 1 S 



27r8 



Ma 



3 + d Vd 
2 + d^¥ 



(— ) 



1+1 



n 



1+1 



(82) 



For comparison with cosmological constraints, it is convenient to determine the graviton 
energy density, normalized to the background photon number density = 2nQ, at the time 
of BBN: 



rrinnQn 



BBN 



TT 



BBN 



45C(3) 



(83) 



The total graviton energy density is then determined by summing over n. For T significantly 
larger than rnxK, we may take the continuum limit 



(84) 



where 



A. 



r 



2, 27r, 47r, 27^^ 



(85) 



for d = 1, 2, 3, 4, . . ., is the area of a unit sphere in d dimensions, and the factor of 1/2"^ in 
Eq. ()84p again accounts for the restriction to n with non-negative components. Integrating 
over all KK levels up to mn = Trh, the total energy density in KK gravitons is 



dn 



V5C(_3) 

271^ 



' /7^,BBN 3 + ^ l/dAg / iRH \ 



2+^ 



(86) 



If KK gravitons decay to SM particles, Cg gives the total energy density, normalized to the 
background photon density, deposited during the era of graviton decay. 

Alternatively, one may take the opposite limit, and assume that KK gravitons decay 
predominantly through KK number preserving interactions —>■ nLKP + X, and nearly 
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all of the energy stored in KK gravitons exists now in the form of KK dark matter. In this 
scenario, the current dark matter energy density, normalized to the critical density, is 



G 



(g 
13 eV 



(87) 



where ~ 410 cm~ is the present photon number density, and ~ 5300 eV cm~ is the 
critical density. 

Numerically, given g^^^ 



3.36 and 0:3 ~ 0.1, we find 



D 



D = 6 



D = 7 



Cg 
^g 
Cg 

^G 

Cg 



1.1 X 10"^^ GeV X C 



200 



8.4 X 10"^ X C 



.200 

8.9 X 10"^^ GeV x C 



1 TeV 



1 TeV 
2 



RH 



RH 



gj^ 
200 



2 



rriKK 

7 
2 



6.8 X 10"^ X C 



.200 

5.0 X 10"^^ GeV X C 



y»KK 
1 TeV 



1 TeV 
1 2 



T, 



RH 



T, 



rrtKK 
n 5 



RH 



200 



mKK 
2 



G 



3.9 X 10^^ X C 



200 



1 TeV 



1 TeV 

2 



RH 



RH 



mKK 



mKK 

2 



(88) 
(89) 
(90) 
(91) 
(92) 
(93) 



Note that for each these expressions, derived from Eqs. (jHUj) and ()87|) . are not mutually 
consistent — the energy density in gravitons is either deposited in SM particles, or in KK 
dark matter (WIMP or superWIMP), or a mixture. These expressions are the maximal 
energy deposited in SM particles, and the maximal primordial energy density of KK dark 
matter. Which, if either, is realized depends on the specific UED model, and detailed 
considerations of higher KK mode cascade decays. On the other hand, these expressions are 
useful, as, if both are within existing bounds, the model is guaranteed to be consistent with 
current constraints. 



VI. COSMOLOGICAL CONSTRAINTS 

There are two sets of constraints on the KK graviton abundance. First, gravitons G" 
may decay to stable LKPs, such as or B^. The fraction of initial energy that winds up 
in LKPs is highly dependent on the UED model and the spectrum of KK modes, which 
determines the form of cascade decays. Clearly the energy density in LKPs cannot exceed 
the energy density of the initial KK gravitons. We may therefore impose the constraint 

^G < ^BM ~ 0.23 . (94) 

If this constraint is satisfied, KK dark matter will not overdose the universe, irrespective of 
details of the UED model. 
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The second set of constraints follows from requiring that SM particles produced in late 
decays of KK gravitons not destroy the successful predictions of BBN. These constraints are 
complicated, depending sensitively on the kind of energy deposited and the time at which it 
is released. For electromagnetic decays, the constraint on = B^^(g has been studied in 
detail 3^ 37, 3^, most recently in Ref. jll]. These constraints are very weak for decay times 
r < 10^ s, increase in stringency from < 10"^ GeV at r ~ 10^ s to C™ ^ lO'^^ GeV at 
r ~ 10^ s, and then remain roughly constant up to r ~ 10^^ s. 

For hadronic cascades, the picture is at present less clear. For r ^ 10^ s, hadronic 
constraints are relatively weak, but they become (q^*^ = B^^'^^g ~ 10^^^ GeV for 10^ s ^ 
r ^ 10^ s 113]. For longer lifetimes, there are no detailed recent analyses. The constraint 
may become weaker, but we assume conservatively that it remains at the (q^'^ ^ 10~^^ GeV 
level. 

From Eqs. (j88|) - (j93p . we see that both the overclosure and BBN constraints may be 
satisfied for any D and Trh ~ ""^kk ~ 1 TeV, even for C ~ 1 and taking the extreme cases 
in which all energy is deposited in either stable LKPs, EM cascades, or hadronic cascades. 
Although there are many uncertainties in this analysis, there is surely an allowed window 
in which the KK graviton abundance satisfies all existing constraints. This is necessary to 
establish the viability of KK WIMP and superWIMP dark matter. For LKPs to achieve 
the desired thermal relic density, either in the form of WIMPs such as or superWIMPs 
such as G^, the universe must be reheated to a temperature Trh ~ mKK/25. The region 
with ideal KK WIMP thermal abundance 0] is also shown. This result shows that there 
are consistent cosmologies in which the reheat temperature is low enough to suppress the 
primordial KK graviton abundance appropriately, but high enough to generate the desired 
thermal relic abundance for WIMP or superWIMP dark matter. 

At the same time, given the rather strong power law dependence of the graviton energy 
density, overclosure and BBN stringently constraint Trh- In Fig. [H we plot the bounds 
on Trh as functions of mKK for D = 5,7 from constraints on overabundance and energy 
release. Roughly we find constraints Trh ^1 — 10 TeV for 100 GeV < mxK < 1 TeV. In a 
given model, one of the constraints may be inapplicable, but they cannot both be completely 
avoided. The range in each constraint results from varying C in the range 0.001 < C < 0.1. 
Note that the Trh limits are rather insensitive to the substantial uncertainties encoded in 
C, given the extreme sensitivity of graviton abundances to Trh; for C varying by two orders 
of magnitude, the bounds on Trh vary by only factors of 2 to 4, depending on the number 
of extra dimensions. 



VII. SUMMARY 

We have provided a general formalism for analyzing the dynamics of gravitons in UED 
theories. In particular, we found the couplings of KK gravitons to fermions and gauge 
bosons and presented the widths for decays of excited fermions and gauge bosons into KK 
gravitons in Eqs. ()61|) and (jMj) . These results are of special relevance when a KK graviton 
is the LKP and a superWIMP candidate, as they determine the observable implications of 
KK graviton dark matter for, for example, BBN, the cosmic microwave background, and 
the diffuse photon flux. 

We have also determined the abundance of KK gravitons produced after reheating in a 
general manner applicable to UED models for arbitrary numbers of extra dimensions, and 
also more generally to other models of extra dimensions. The possibility of populating a large 
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0.1 0.2 0.5 1.0 2.0 5.0 10.0 



rriKK (TeV) 

FIG. 1: Bounds on the reheat temperature Trh as a function of mxK from the overclosure constraint 
< 0.23 and the BBN constraint C,g < 10^^^ GeV for D = 5,7, as indicated. We assume 
g^^ = 200; the range in each bound arises from varying C from 0.001 to 0.1 (see text). In the 
region with Trh < '"T'Kk/25, Trh is too low to generate the thermal relic abundance for WIMPs. 
The vertical bands delimit regions where the thermal relic abundance is too low (r^wiMP < 0.1), 
approximately right, and too high (f^wiMP > 0.3). 



number of graviton states at different KK levels implies that the production of gravitons 
after reheating is extremely efficient and extremely sensitive to the reheat temperature Trh- 
For d extra dimensions, the energy density in gravitons, given in Eq. ()86|) . is 



Cg = 2^ m^n — . 95) 

Y s M4 VmKK/ 

This is to be contrasted with the case of gravitinos [2^ 2^, for which 

^- =^-!^ ^ !!!iusY ^RH_ /gg^) 

^ ^ s M4 msusY ■ ^ ^ 

The constraints on Trh are therefore extremely stringent. They are presented in Fig. and 
are of the order of Trh < 1 - 10 TeV for 100 GeV < tukk < 1 TeV. 

The constraints derived here are robust, being independent of the gravi-scalar mass js^ 
and applicable irrespective of which KK particle is the LKP. These constraints also apply 
in the presence of KK-parity violating interactions, as these will only serve to increase the 
primordial graviton production and lead to the decay of gravitons to SM particles. They 
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supplement the req uirement that the reheat temperature be below the cutoff of the 4D 
effective theory j35f] and rather severely constrain ideas for leptogenesis. Such low reheat 
temperatures also constrain inflation scenarios, requiring, for example, that inflaton decay to 
SM particles be suppressed by extremely small couplings or kinematically through enhanced 



SM plasma masses at high temperatures |41 



At the same time, we have found that there exists a range of reheat temperature with 
Trh ~ "^KK such that the primordial production of KK gravitons is within bounds, but a 
thermal relic density of WIMPs may be produced. KK WIMP and superWIMP candidates 
are therefore viable, despite the stringent graviton constraints applicable to these extra 
dimensional theories. 
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